Chiral three-nucleon forces and bound excited states in neutron-rich oxygen isotopes 
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We study the spectra of neutron-rich oxygen isotopes based on chiral two- and three-nucleon 
interactions. First, we benchmark our many-body approach by comparing ground-state energies 
to coupled-cluster results for the same two-nucleon interaction, with overall good agreement. We 
then calculate bound excited states in 21 > 22 > 23 0, focusing on the role of three-nucleon forces, in 
the standard sd shell and an extended sdf 7 / 2 p3/2 valence space. Chiral three-nucleon forces provide 
important one- and two-body contributions between valence neutrons. We find that both these 
contributions and an extended valence space are necessary to reproduce key signatures of novel shell 
evolution, such as the N = 14 magic number and the low-lying states in 21 and 23 0, which are 
too compressed with two-nucleon interactions only. 

PACS numbers: 21.60.Cs, 21.10.-k, 27.30.+t, 21.30.-x 



Introduction. - The oxygen isotopes provide an excit- 
ing laboratory, both experimentally and theoretically, to 
study the structure of extreme neutron-rich nuclei to- 
wards and beyond the neutron dripline at 24 Since 
24 O is anomalously close to the valley of stability, all 
bound oxygen isotopes have been explored experimen- 
tally ■ The neutron- rich isotopes exhibit surprising 
new features, such as the closed-shell properties of 22 
and 24 O, and the corresponding appearance of new magic 
numbers N = 14 and N = 16 [|, \M^. The spectra with 
these fingerprints provide a challenge for nuclear theory, 
and a unique test of nuclear forces and many-body meth- 
ods for neutron-rich nuclei. 

While the chain of oxygen isotopes can be explored by 
different theoretical methods fl(ML4l ] , recent studies have 
focused on ground-state properties relevant for under- 
standing the location of the neutron dripline. In theories 
based on two-nucleon (NN) forces only, the interactions 
were found to be too attractive and to predict bound 
oxygen isotopes to 28 [lfj. Coupled-cluster (CC) cal- 
culations obtained bound isotopes to 25 O with chiral NN 
forces only, but found the location of the dripline to de- 
pend on the resolution scale, implying that higher-body 
forces need to be considered [ll| . The first work including 



three-nucleon (3N) forces demonstrated that they play a 
decisive role for the oxygen anomaly and can exp lain why 
24 is the heaviest bound oxygen isotope [lOj . Subse- 
quent CC calculations implementing adjusted 3N effects 
and coup ling to the continuum have found similar re- 
sults [II]. 

Calculations of excited states in neutron-rich oxygen 
isotopes, when based on NN forces only, show similar de- 
ficiencies, evident in the poor description of spectra (see 
Figs. [SHU, where low- lying states arc too compressed. 
In 23 0, CC calculations based on NN forces [1J] found 
the first excited | state to be more than 2 MeV below 
experiment. As a result, phenomenological shell-model 
interactions fl5l - [l7l ] have been developed that are ad- 
justed to correct for these deficiencies. As suggested by 



Zuker [18(, we have shown that these phenomenological 
adjustments can be largely traced to neglected 3N forces, 
both for the oxygen ground-state energies [l(| EH and in 
calcium where 3N forces are key for the N = 28 magic 
number (20[. Excited states in the oxygen isotopes have 
also been studied recently neglecting 3N forces [HJ and 
in 18 - 19 with density-dependent 3N forces [23|, but both 
rely on empirical inputs. 

Three-nucleon forces, first introduced in the pioneering 
work of Fujita and Miyazawa [24| , arise because nucleons 
are finite-mass composite particles that can be excited 
by interacting with other particles. The long-range part 
of 3N forces is dominated by the Fujita and Miyazawa 
singlc-Dclta-cxcitation mechanism, where one nucleon 
virtually excites a second nucleon to the A(1232McV) 
resonance that is then de-excited by interacting with a 
third nucleon. Additional long- and shorter-range 3N in- 
teractions are included naturally in chiral effective field 
theory (EFT) (25|, which provides a systematic expan- 
sion for nuclear forces. Similar to the oxygen and calcium 
isotopes, chiral 3N forces give repulsive contributions to 
the neutron and nuclear matter energy [26T - |28| and to 
pairing gaps in nuclei (29| . 

In this paper, we present the first study of the impact 
of chiral 3N forces on bound excited states in neutron-rich 
oxygen isotopes. We first discuss the forces and many- 
body approach used, and benchmark the ground-state 
energies with CC theory based on the same NN forces. 
We then calculate the spectra of 21 ' 22 > 23 in the stan- 
dard sd shell and an extended sdf7/ 2 p3/ 2 valence space, 
focusing on the role of 3N forces. Our results show that 
one- and two-body contributions from chiral 3N forces to 
valence neutrons, included to third order in our many- 
body formalism, as well as a valence space larger than 
the standard sd shell are important to understand the 
experimental excited state structures. 

Interactions among valence neutrons.- We calculate 
the interactions among valence neutrons based on chi- 
ral NN and 3N forces following Refs. @, [H|. At the 
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FIG. 1. Ground-state e nerg ies of neutron-rich oxygen isotopes relative to O, with experimental energies from the AME 2003 
atomic mass evaluation [3J]. Panel (a) compares results calculated in MBPT to those from CC theory based on the same low- 
momentum NN interactions. The MBPT results are based on the SPEs obtained in CC theory, for details see text. Panel (b) 
shows the NN-only MBPT results at 2nd and 3rd order in the sd shell, and 3rd order in the sdf 7 / 2 p3/ 2 shell, using empirical 
SPEs. In panel (c), we present the energies obtained at 3rd order from NN and 3N forces, with calculated MBPT SPEs of 
Table [I] that also include 3N contributions, in comparison to the phenomenological USDb model. 



NN level, we first perform a renormalization-group (RG) 
evolution [13 of the N 3 LO NN potential of Ref. [3l| to 
a lower resolution scale A = 2.0 fm . The resulting 
low-momentum interaction Vj ow k is used to calculate the 
two-body interactions among valence neutrons to third 
order in many-body perturbation theory (MBPT), fol- 
lowing the formalism of Ref. (32[. We use a harmonic 
oscillator basis of 13 major shells with huj = 13.53 McV, 
appropriate for the oxygen isotopes. In the MBPT cal- 
culations including 3N forces, we also scale the resulting 
valence-shell interaction and bound single-particle ener- 
gies (SPEs) by hjj ~ A" 1 / 3 . The MBPT results are 
converged in terms of intermediate-state excitations. In 
addition to the standard sd shell, we also consider an ex- 
tended sdf 7 /jp 3 / 2 valence space, where the largest changes 
are due to including the lower-lying / 7 / 2 orbital. 

To study the validity of the MBPT approach, we have 
carried out CC calculations for the ground-state ener- 
gies of the oxygen isotopes, using the same Vi ow fc in- 
teraction and the same basis space of 13 major shells 
with ftuj = 14 McV. The CC energies are well con- 
verged, with a flat hut dependence of ~ 200 keV change 
from hw = 14 — 20MeV. The results are shown in 
Fig.[T](a) relative to the ground-state energy of 16 0. The 
closed j'-subshell systems, 16 ! 22 i 24 > 28 o, are calculated at 
the A-CCSD(T) level (with small differences from 
CCSD ~ lMeV). The A± 1 systems, i7,2i,23,25,27 0) arc 
obtained with the CC particle- attached/removed equa- 
tions of motion method at the singles and doubles level 
(PA/PR-EOM-CCSD) [H|. 

For the CC comparison, wc perform the MBPT cal- 
culations in the sd shell, where for consistency the SPEs 
are taken as the PA-EOM-CCSD (d 5 / 2 , su 2 , d 3 / 2 ) energies 
in 17 0. The particle-attached fy 2 and p a / 2 energies are 
not single particle in character (at this NN level), so the 
MBPT s<i-shcll comparison provides the cleanest bench- 



mark. In Fig. [T] (a), we find that the MBPT energies 
agree within a few percent with CC theory (the agree- 
ment of excitation energies is expected to be even better). 
This shows that, at this level, MBPT can be compara- 
ble to CC theory for Vj ow fc interactions, but also high- 
lights the important role of SPEs. Furthermore, compar- 
ing with the results shown in Figs. Q] (b) or (c), we see 
that the addition of 3N forces or going to an extended 
valence space has a greater impact on the ground-state 
energies than the uncertainty in MBPT in Fig. [1] (a). 

At the 3N level, we include 3N forces among two va- 
lence neutrons and one nucleon in the core. This corre- 
sponds to the normal-ordered two-body part of 3N forces, 
which was found to dominate in CC calculations [35j 
over residual 3N forces (the latter contributions should 
be weaker for normal Fermi systems I36|). We take 
into account chiral 3N forces at N 2 LO (37| that give 
rise to repulsive monopole interactions between valence 
neutrons [13, H3|. At N 2 LO, this includes long-range 
two-pion-exchange parts a (due to A and other exci- 
tations), plus shorter-range one-pion exchange cu and 
3N contact eg interactions, with the cd,ce couplings 
fit to the 3 H binding energy and the 4 He radius for 
A = A 3N = 2.0 fm" 1 [Hi. The fit to light nuclei ap- 
proximately includes effects of the RG evolution in the 
3N sector, up to higher-order many-body forces. For all 
results, these 3N forces, within 5 major oscillator shells, 
are included fully to third order in MBPT. 

For the SPEs in 17 0, we solve the Dyson equation, con- 
sistently including one-body contributions to third order 
in MBPT in the same space as the two-body interactions. 
While the convergence in terms of intermediate-state ex- 
citations is slower for the one-body energies, the results 
are still well converged in 13 major shells. In the Dyson 
equation, we include chiral 3N forces between one va- 
lence neutron and two core nucleons, which corresponds 
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TABLE I. Empirical and calculated (MBPT) SPEs in MeV. 



orbital 


sd shell 
USDb [15] MBPT 


sdf 7/2 p 3/2 


shell 
MBPT 


^5/2 


-3.93 


-3.78 


-3.95 


-3.46 


Sl/2 


-3.21 


-2.42 


-3.16 


-2.20 


da/a 


2.11 


1.45 


1.65 


1.92 


jV/2 






3.10 


3.71 


P3/2 






3.10 


7.72 



to the normal-ordered one-body part of 3N forces. The 
3N contributions range between 1.5 — 3.0 MeV, approxi- 
mately an order of magnitude larger than typical normal- 
ordered two-body matrix elements. This is consistent 
with the normal-ordering hierarchy (35l . [36j . At each 
self-consistency iteration, the unperturbed harmonic os- 
cillator spectrum is updated such that the energy of the 
degenerate valence-space orbitals is set to the centroid, 
Ej(2j + + 1), of th e calculated SPEs of the pre- 

vious iteration. 

Our results for the SPEs in the sd and sd]i/ 2 p3/ 2 va- 
lence spaces are given in Table HI in comparison with 
empirical SPEs taken from the USDb QJ] and SDPF- 
M [16| models. With 3N forces included, the resulting 
SPEs are similar to the empirical values, except that the 
p3/ 2 orbit is significantly higher. Because the higher-lying 
orbitals are unbound, a necessary future improvement 
will be to represent these correctly as continuum states. 
Work in this direction is under development, but CC cal- 
culations in this region have found that continuum cou- 
pling effectively lowers energies of unbound states only 
by - 300 keV 

We finally investigate possible center-of-mass (cm) 
contamination by adding a cm term, /3H cm , with /3 = 0.5, 
to our original Hamiltonian. While the qualitative pic- 
ture remains, there is some additional compression in the 
studied spectra as a result. This is not unexpected since 
the non-zero cm two-body matrix elements are also rele- 
vant matrix elements of the MBPT calculation, and mak- 
ing a clear separation between them is therefore difficult. 
A better test of cm factorization has recently been pro- 
posed in the context of CC theory [39| , which we plan to 
adapt to the shell model context in the future. 

Ground-state energies.- Although the focus of this 
work is on excited states, we first discuss results for the 
ground-state energies in this framework. All calculations 
based on NN forces are performed with empirical SPEs, 
and all calculations including 3N forces with the consis- 
tently calculated MBPT SPEs of Table Q] We take into 
account many-body correlations by diagonalization in the 
respective valence spaces. Figures [T] (b) and (c) show 
the resulting ground-state energies compared with ex- 
periment. With NN forces only, all isotopes out to 28 O 
remain bound because the two-body interactions among 
valence neutrons are too attractive [Toj j . Going to the 
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FIG. 2. Excitation energies of bound excited states in 1 
compared with experiment 0, The NN-only results are 
calculated in the sd and sdfr/ 2 p3/ 2 shells with empirical SPEs. 
The NN+3N energies are obtained in the same spaces, but 
with calculated MBPT SPEs, including 3N contributions to 
valence two-body interactions and SPEs. The dashed lines 
give the one-neutron separation energy S n . 



extended sdf r / 2 py 2 space provides additional attraction, 
mainly due to the correlations involving the f 7 / 2 orbital. 
To gauge the order-by-order convergence of MBPT, we 
also show in Fig. [1] (b) the ground-state energies to 2nd 
order. The small difference from 2nd to 3rd order is con- 
sistent with the CC agreement and shows that the MBPT 
calculations can be converged to a good level. 

When 3N forces are included in Fig.[T](c), we see their 
clear repulsive contribution to the ground-state energies, 
relative to 16 In both valence spaces, there is sig- 
nificantly less binding compared to the NN-only calcula- 
tions. The calculated energies in the sd shell are under- 
bound with respect to experiment, while in the extended 
sdf r / 2 p3/ 2 space, we find very good agreement. Finally, 
it is interesting to note that in the sdfr/ 2 p3/ 2 space, the 
ground state wavefunctions contain a significant admix- 
ture of a 2p2h intruder configuration. The occupancies 
of the ground states of 22 O and 24 O are as expected from 
the sd-shell except for a prominent fy 2 component al- 
ready in these isotopes. The fi/ 2 occupancies are 1.3 and 
1.6 for 22 and 24 0, respectively (with the ds/ 2 occupan- 
cies reduced to 4.2 and 4.4). This configuration largely 
persists in the excited states discussed below. 

Spectra.- Next, we study the properties of bound ex- 
cited states in the neutron-rich oxygen isotopes 21 > 22 > 23 0, 
which have been the subject of recent experimental in- 
vestigations @, 0, H, [|| . While we will not explicitly dis- 
cuss spectra of the lighter isotopes, we note that in the 
NN+3N extended-space calculation, the first excited 2 + 
(4+) energies in 18 and 20 O arc 2.25 MeV (3.90 MeV) 
and 2.14 MeV (4.99 MeV), respectively, compared with 
experimental values 1.98 MeV (3.55 MeV) and 1.67MeV 
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FIG. 3. Excitation energies of bound excited states in 22 O 
compared with experiment ||, U (labels are as in Fig. [2]| . 



FIG. 4. Excitation energies of bound excited states in 23 
compared with experiment [f| (labels are as in Fig. [2}. 



(3.57 MeV). Furthermore, the calculated E(4+)/E(2 + ) 
ratios of 1.79 and 2.33 for 18 and 20 O compare well 
with the experimental values of 1.73 and 2.14. For each 
nucleus we discuss and compare with all measured states 
below the one-neutron separation energy S n - The cou- 
pling to the scattering continuum is important for un- 
bound states and as the energies of excited states ap- 
proach s n [HE 53]. 

Here, we focus on 3N force con- 
tributions and note possible continuum effects for the 
higher-lying states. 

21 0.~ The spectrum of 21 in Fig. [5] shows clearly the 
deficiencies of calculations based on NN forces. In a basic 
sd-shell picture, the ground state is interpreted as a one- 
hole | + state below the 22 O closure, while the excited 
states result from a one-particle excitation to the Si/ 2 
orbital coupled to two holes in d$/ 2 with angular momenta 

+ , 2 + , and 4 + (a | + state is observed above the neutron- 
decay threshold). These states therefore probe the d 5 / 2 — 
Si/ 2 gap. With NN forces only, the calculated spectrum 
is in very poor agreement with experiment, being too 
compressed for the low-lying states and too separated 
for the | + and | + states. Expanding the valence space 
to the sdfii 2 p3/ 2 shell provides only slight improvement 

for the | state, but otherwise has little effect. 

With the inclusion of 3N force contributions, there is 
a clear improvement in spacing already in the sd shell, 
although with an incorrect ^ + ground state. For the sd 
shell, S n is well below the experimental value, which re- 
flects the undcrbound ground-state energies. Our results 
in the extended sdf7/ 2 p3/ 2 space with 3N forces show sig- 

1 + q + 

nificant improvements, although the ^ and ^ states re- 
main somewhat below the experimental values. Further- 
more, due to their presence above the continuum thresh- 
er _)_-_(_ 

old, the | 2 and ^ energies may be reduced, bringing 
them closer to experiment. 



22 0.- The first oxygen isotope exhibiting closed shell 
properties for a non-standard magic number, O, has its 
first 2 + state at almost twice the energy as those in 18 
and 20 O. In contrast to 24 0, whose closed-shell nature 
can be qualitatively well described in NN-only calcula- 
tions due to the large separation between the g?3/ 2 and 
Si/ 2 SPEs (see, e.g., Ref. [Hj]), the spectrum of 22 is not 
well reproduced with NN forces, as shown in Fig. [3] In 
the scf-shell, the ground state is predominantly a filled 
d&/ 2 subshell, while the first excited 2+ and 3 + states re- 
sult from a one-particle Si/ 2 excitation coupled to a ds/ 2 
hole. The higher- lying + , 2 + , 4 + triplet is formed by 
two excited su 2 particles coupling to two ds/ 2 holes. 

With NN forces in the sd shell, the first 2 + state is 
below experiment, the spacing to the 3 + state too broad, 
and the + , 2 + , and 4 + states are too compressed. Ex- 
tending the space provides only a slight improvement 
compared to the experimental spectrum. The addition 
of 3N forces in the sd shell reduces the first 2 + energy 
but otherwise improves the level spacings. 

It is only when we move to the sdfy 2 py 2 space with 
3N forces included that we see significant improvement 
in the spectrum, with the 2+ energy now being approx- 
imately twice the 2+ energies of 18 and 20 O. Although 
the 2+ - 3+ splitting is reduced, it is still somewhat large 
and results in an inversion of the 3 + and + states com- 
pared with experiment, a feature also seen in all phe- 
nomenological USD models. Since the 4 + energy lies 
above S n , some lowering can be expected from continuum 
effects. Similar to N — 28 in calcium [2(j, we find that 
both 3N forces and an extended valence space are im- 
portant for reproducing properties of the N = 14 magic 
number in oxygen. 

23 O.- Occupying a position of particular interest be- 
tween two doubly-magic isotopes, the spectrum of 23 
provides a unique test for theory, as it should simulta- 
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ncously reflect the features of both. The se?-shell 23 O 
ground state is dominated by one particle in the sy 2 or- 
bit, while the two lowest excited states are expected to 
be a single particle | + one- hole excitation, indicative of 
the strength of the 22 O shell closure, and a higher- lying 
single particle ^ one-particle excitation, reflecting the 
strength of the 24 shell closure. Recent experimental 
measurements of the lowest-lying states in 23 O shown in 
Fig. |4] confirm this picture, although no bound excited 
states are found |f|. The | + state lies just above the 

neutron decay threshold, while the ^ state resides well 
in the continuum, and a complete description cannot ig- 
nore these effects. 

From Fig. [4] we see that the | state is bound in all 

calculations. With NN forces in the sd shell, the | + 

i + ^ + 

state even lies below the ^ , and the % is bound, con- 
sistent with the expectation of a bound dy 2 orbital at 
24 without 3N forces Extending the space to the 
sdfr^ 2 p 3 / 2 orbitals is a clear improvement, and the | is 

now 0.54 MeV above the \ + ground state. This is similar 
to the 0.35 MeV excitation energy obtained in CC theory 
based on a N 3 LO NN potential only 

The repulsive 3N contributions improve the spectrum 
considerably. We again find the best agreement with ex- 
periment in the extended sdfi/ 2 p 3 / 2 space, where the en- 
hanced d5/ 2 — fr/2 correlations, important for the high 2 + 

in 22 O, bring the | + energy very close to the experimen- 
tal value. This shows that 3N forces are essential for the 
| state, increasing the excitation energy by ~ 2 MeV 

compared to the NN-only result. The | state, however, 
is significantly lower than experiment lying just beyond 
the neutron-decay threshold. When continuum coupling 



is included, this will effectively lower the d 3 / 2 orbit 
pointing to the need to explore uncertainties as well as 
further improvements in our theoretical approach. 

Summary.- We have presented the first study of the 
impact of chiral 3N forces on the spectra of neutron-rich 
oxygen isotopes. Based on the MBPT approach and chi- 
ral EFT interactions, we find that the traditional sd shell 
is insufficient, even with 3N forces, to account for key 
experimental features in the spectra, in particular the 
N = 14 magic number. Going to an extended sdf7/ 2 p 3 / 2 
valence space, and including NN and 3N contributions 
to third order in MBPT consistently to the two-body 
interactions and the one-body SPEs, our results are in 
good overall agreement with experiment. It will be in- 
teresting to study how the physics of the extended va- 
lence space may be renormalized into phcnomenological 
s<i-shell models. In addition, extending our framework to 
include continuum effects for weakly-bound or unbound 
states, and for nuclei beyond the driplines, as well as de- 
veloping nonperturbative approaches to valence shell in- 
teractions (see, e.g., Refs. [Lu443l ]) are important future 
steps. 
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